Measuring Double Parton Distributions in Nucleons at 
Proton-Nucleus Colliders 



M. Strikman 

Department of Physics, Pennsylvania State University, University Park, PA 16802 

O ' D.Treleani 
O 

' Universitd di Trieste; Dipartimento di Fisica Teorica, Strada Costiera 11, Miramare-Grignano, 

> : 

, ^ ■ and INFN, Sezione di Trieste, 1-34014 Trieste, Italy. 

o 

cn 

^ , Abstract 
> 
oo 

Tlj- ' We predict a strong enhancement of multijet production in proton-nucleus 

collisions at collider energies, as compared to a naive expectation of a cross 
section oc A. The study of the process would allow to measure, for the first 
time, the double parton distribution functions in a nucleon in a model in- 
dependent way and hence to study both the longitudinal and the transverse 
correlations of partons. 
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The systematic studies of hard inclusive processes during the last two decades have let 
to a pretty good understanding of the single parton densities in nucleons. However very 
little is known about multiparton correlations in nucleons which can provide a new tool for 
discriminating between different models of nucleons. Such correlations may be generated, for 
example, by the fluctuations of the transverse size of the color field in the nucleon leading, 
via color screening, to correlated fluctuations of the densities of gluons and quarks. A related 
source of correlations is QCD evolution, since a selection of a parton with a given x, may 
lead to a local (in transverse plane) enhancement of the parton density at different x values. 

It was recognized already more than two decades ago [|^ that the increase of parton 
densities at small x leads to a strong increase of the probability of nucleon-nucleon colli- 
sions where two or more partons of each projectile experience pair-vice independent hard 
interactions. Although the production of multijets through the double parton scattering 
mechanism was investigated in several experiments at pp,pp colliders, the interpreta- 
tion of the data was hampered by the need to model both the longitudinal and the transverse 
partonic correlations at the same time. The aim of this letter is to point out that the near 
future perspectives to study proton-nucleus collisions at RHIC, as well as the plans for 
pA collisions at LHC, provide a feasible opportunity to study separately the longitudinal 
and transverse partonic correlations in the nucleon as well as to check the validity of the 
underlying picture of multiple collisions. 

The simplest case of a multiparton process is the double parton collision. Since the 
momentum scale pt of a hard interaction corresponds to much smaller transverse distances 
~ 1/pt in coordinate space than the hadronic radius, in a double parton collision the two 
interaction regions are well separated in the transverse space. Also in the cm. frame pairs of 
partons from the colliding hadrons are located in pancakes of thickness < (l/a;i + l/x2) /Pc.m.- 
So two hard collisions occur practically simultaneously as soon as Xi^X2 are not too small and 
there is no cross talk between two hard collisions. A consequence is that the different parton 
processes add incoherently in the cross section. The double parton scattering cross section, 
being proportional to the square of the elementary parton-parton cross section, is therefore 
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characterized by a scale factor with dimension of the inverse of a length squared. The 
dimensional quantity is provided by the nonperturbative input to the process, namely by the 
multiparton distributions. In fact, because of the localization of the interactions in transverse 
space, the two pairs of colliding partons are aligned, in such a way that the transverse 
distance between the interacting partons of the target hadron is practically the same as the 
transverse distance between the partons of the projectile. The double parton distribution 
is therefore a function of two momentum fractions and of their transverse distance, and it 
can be written b). Actually F depends also on the virtualities of the partons, 

Q^,Q'^, though to make the expressions more compact we will not write explicitly this 
dependence. Hence the double parton scattering cross section for the two "two two" 
parton processes a and (3 in an inelastic interaction between hadrons a and b can be written 
as: 

a£){<y,P) — ^ J ^ai^i, X2]b)aaXi, x[)af3{x2, x'2)Tf,{x[, X2;b)dxidx[dx2dx'2(fb (1) 

where m = 1 for indistinguishable parton processes and m = 2 for distinguishable parton 
processes. Note that, though the factorization approximation of Eq.(|l|) is generally accepted 
in the analyses of the multijet processes and appears natural based on the geometry of the 
process, no formal proof exists in the literature. As we will show below the study of the 
A-dependence of this process will allow to perform a stringent test of this approximation. 

In the case of A^A^ scattering one cannot proceed further without making some simplifying 
assumptions about transverse correlations of partons in nucleons. Our key observation is that 
the introduction of a new large transverse scale: the nucleus radius, allows to separate the 
effects of the transverse and longitudinal parton correlations. Essentially, we can express the 
function TA{xi,X2,b) through Tiy{xi,X2,b) and the distribution of nucleons in the nucleus, 
practically without any extra model assumption. Here to simplify the discussion we neglect 
small non-additive effects in the parton densities, which is a reasonable approximation for 
0.02 < X < 0.5. In this case we have to take into account only b- space correlations of 
partons in individual nucleons. One has therefore two different contributions to the double 
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parton scattering cross section. The first one, af , which is represented in Fig, la, is the same 
as for the nucleon target (the only difference being the enhancement of the parton flux) and 
the corresponding cross section is 

af = f^D y <fBT{B) = Aan. (2) 

where T{B) is the nuclear thickness, as a function of the impact parameter of the hadron- 
nucleus collision B. 

The contribution to the term in Ga {x'l ,X2,b) due to the partons originated from different 
nucleons of the target (Fig.l.b) can be calculated solely from the geometry of the problem 
by observing that the nuclear density does not change within a transverse scale (b) ^ Ra- 

The two simplest methods are to use the AGK cutting rules 0, or the technique of 0. 
We can write, for two indistinguishable parton processes, 
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G]\f{xi,X2)a'{xi,x[)a{x2,X2)GN{x[)Gi\i{x2)dxidx[dx2dx2 / d BT (B), (3) 



where Gn{xi,X2) = J d'^bTN{xi,X2;b), while are nucleon and x'^ are nuclear parton frac- 
tions. The distinctive feature of the ct2- term is in difference from the case of A^A^ interactions, 
since no transverse scale factor related to the nucleon scale is present in cr|^. The correct 
dimensionality is provided by the nuclear thickness function, which appears in a2 at the 
second power. The two contributions and o"!' are therefore characterized by a different 
dependence on the atomic mass number of the target. The A-dependence of the two terms 
is in general a function of the values of the momentum fractions and of the virtuality scale of 
the 2 — s> 2 interactions. The simplest situation is in the kinematical regime where shadowing 
corrections to the nuclear structure function can be neglected, erf is then proportional to 
and af to A^'^. (Note that the nuclear surface effects lead to a faster dependence of 
/ T'^{B)d'^B on A, for A < 240, than the naive expectation A^/^). The presence of two terms 
with distinctive A-dependence (and comparable magnitude for a wide range of x, see Eq.(|]) 
below) will allow to separate them easely experimentally and also to check in the course of 
such an analysis the factorization approximation of Eq.(|l]). 



To estimate the relative importance of erf and cr|^, and only to that purpose, we 
use the CDF analysis [Q] were it was assumed that all correlations in fractional mo- 
menta may be neglected, so that one can write Fat in a factorized form as a product of 
two parton densities Gn{x) and of a function of the inter-parton transverse distance b: 
Tn{x,x' ,b) = G]\f{x)Gj\f{x')F{b). (In a sense this could be considered as marely a conve- 
nient parametrization of the experimental data). With these simplifications one obtains 

aD{a,l3) = — , (4) 

where cTq, and are the inclusive cross sections for the two processes a and f3 in the hadron 
- hadron interactions. Under the factorization assumption the whole new information on 
the hadron structure can be reduced to a single quantity with dimensions of a cross section, 
(Te// which was measured by CDF to be: 

(Te// = 14.5± 1.7+^:5m6. (5) 

Within the accuracy and in the limited kinematic range accessible to the experiment 
(0, 01 — 0.40 for the photon+jet scattering, 0.002 — 0.20 for the dijet scattering) no evidence 
was found of a x dependence of cie//, supporting the simplest uncorrelated picture of the 
interaction. However the absolute value of the cross section is significantly larger (by a 
factor > 2) than the cross section one would obtain within the factorization hypothesis, 
when assuming that the transverse distribution of partons reflects the matter distribution 
in the nucleon needed to obtain the value of the nucleon non-single-diffractive cross section 
measured by CDF (for extensive discussions see @] ). So the CDF data actually appear 
to indicate the presence of parton-parton correlations in a nucleon, though one cannot 
distinguish whether they are solely due to transverse correlations or to a combination of 
longitudinal and transverse correlations. 

Substituting Eq.(4) in Eq.(2) and the factorized form of Fat in Eq.(3) we can estimate 
the relative importance of erf and af : 

J, = ^J^-=//«0.45. -) (6) 
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Here we evaluated / T'^{B)(PB by using the standard experimentally determined Fermi 
step parameterization of the nuclear matter densities ||^. One can see from Eq.^ that for 
heavy nuclei, which are available at RHIC, the second term will constitute about 70% of 
the cross section and hence the study of the A-dependence of the four jet production will 
allow a straightforward separation of the two contributions to the cross section. (Note 
also that the cross section of the two partons four jets process, which constitutes a 
background to the four— four processes, depends linearly on A, so that its contribution may 
be disentangled by studying the A-dependence of the cross section). It is worth emphasizing 
that, if the small value of (Je// is due to the correlation of the longitudinal distributions, the 
relative contribution of the second term would be further enhanced. Considering that the 
enhancement of o"/) in pp collisions is roughly by a factor two as compared with the naive 
expectation, one would expect in this case an additional enhancement of cr^ by a factor 
~ ^2. 

One could question whether the soft particle production background may create more 
serious problems in pA scattering than in A^A^ scattering. It appears that this problem can 
be avoided by choosing Xi > x'^, Xi > 0.1 and selecting a kinematics close to 90 degrees in 
the cm. of the partonic collisions. In this case the jets are produced predominantly in the 
proton fragmentation region, where the soft hadron multiplicity in pA collisions is smaller 
than in A^A^ collisions. 

In order to extend the analysis to the x'^ < 0.01 kinematics one needs to take into 
account shadowing effects in the nuclear parton densities. Here we restrict our discussion 
to the case of the leading twist parton shadowing, which is a pretty safe approximation 
for pt > 5 — 7GeV/c and x > 10^^, which is anyway the minimal cut on pt to be able to 
observe the jets. (For recent estimates of the kinematics were black body/unitarity effects 
may become important see |]^). In the case where only one of the nuclear partons is in the 
shadowing region, the ratio of cr|^/o"f is modified only by the dependence of the shadowing 
on the nuclear impact factor, given the different selection of impact parameters in two terms 
af and o"^ (the integration with measure T'^{B) leads in fact to a somewhat smaller average 
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B as compared with the integration with measure T{B)). We have performed a numerical 
estimate of this effect within the leading twist approximation of ^ and we found that even 
for A = 240 this effect leads to a decrease of / erf by < 10%. The effect could be in any case 
studied experimentally by investigating the single hard scattering as a function of A in the 
same kinematics. When both x[s are in the shadowing region the evaluation of the effect is 
more model dependent, though it still appears to be rather small. In any case such kinematics 
is more appropriate for the study of the dynamics of the nuclear shadowing and hence is not 
directly related to the subject of this letter. Note also that pushing such measurements into 
a kinematical region close to the black body limit would create additional problems since, 
due to the increase of the transverse momenta of the nuclear partons, the pair- vice azimuthal 
correlation, which allows an easy identification of the double parton collision events, would 
become weaker and weaker. Note also that we argued in the beginning that processes which 
may violate factorization should be amplified when Xi,x'^ become smaller. Hence the study 
of the A-dependence in this kinematics would provide an additional test of the factorization 
approximation. 

Summarizing, the study of the A-dependence of the double parton scattering will allow 
to separate two contributions to the cross section - due to scattering off one and two nucleons 
of the nucleus. Because of the large nuclear size: Ra ^ Rn, the erf term provides a model 
independent measurement of the double parton densities in nucleons - while no such model 
independent measurement is possible with proton targets. At the same time the comparison 
of the two terms will allow a practically model independent determination of the transverse 
separation between two partons (modulus a possible small effect due to a different transverse 
separation of r(xi, X2, b) and r{x'i,X2, b)) as well as checking the factorization approximation. 

Obviously one can consider also three parton collisions. In difference from the case of 
the double collisions it is more difficult for the available range on nuclei to extract the 
triple parton distribution without making simplifying hypotheses. This is because the triple 
scattering process originates due to three different mechanisms, corresponding to the number 
of target nucleons involved. While the terms with one (Fig. 2a) and three target nucleons 
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(Fig. 2c) are analogous to the contributions already considered for the double scattering, 
the contribution with two different target nucleons - Fig. 2b, is different. In the latter case, 
in fact, the integration on the transverse coordinates of the interacting partons involves at 
the same time two partons of the projectile and two partons of the target. The simplest 
possibility is that the longitudinal and transverse degrees of freedom can still be factorized, 
in this case the integration over the transverse partonic coordinates gives as a result a factor 
with dimensions of the inverse of a cross section. We call the new dimensional quantity 
a'gff and a naive expectation would be that its value is not much different from (Jeff- The 
different contributions to the triple scattering cross section are therefore: 



where ut is the triple parton scattering cross section on a nucleon target. The second 
term provides an additional information about correlations of partons in nucleons while the 
third term measures triple parton density in nucleons. If we assume that the integral over 
transverse coordinates in dimension scale for erf, a J are approximately and a~jj: we can 
estimate that the relative importance of the three terms, for A > 10, is approximately: 



This estimate indicates that the A-dependence of is much stronger than for a^, with 
the scattering off several nucleons becoming important already for light nuclei. The crj 
term is likely to become comparable to the other terms for heavy nuclei, so in principle 
an accurate study of the A-dependence would allow to measure all three terms separately 






1 : 1.45(^/10)°-^ : 0.25(^/10) 



(8) 
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and hence determine the triple parton density in a nucleon in a model independent way. 
Obviously one would need LHC energies and a large acceptance FELIX-type detector 
to be able to study such reactions. 

One can go a step further and try to get information about global characteristics of the 
nucleon as a function of the values of the flavor, x's, etc of the probed partons. Really 
the number of the secondaries produced gives indications of the actual transverse size of 
the projectile, so that one would expect to observe a relatively smaller population of sea 
quarks and gluons in events with few secondaries, while, in that case, the momentum carried 
by the valence should be larger than average. Hence one might start by studying the 
correlation between soft characteristics of the events (the simplest for RHIC would be number 
of neutrons in the zero angle calorimeter) and the momentum fraction x of the projectile 
parton in a single hard collision |jlO[. A next step would be to compare the single and 
double hard scattering events for fixed Xi, say xi ~ 0.2, while increasing the value of X2- 
If, when selecting two fast partons in a nucleon, one selects configurations with a small 
size, one would expect, for example, that the number of knockout nucleons would decrease 
when X2 increases. Overall, such studies would allow to obtain unique information about 
the three-dimensional structure of the nucleon. 
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FIG. 1. Two contributions to the "four four" process in pA scattering. Dashed Hues represent 
hard interactions. 
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(a) (b) (c) 

FIG. 2. Three contributions to the "six six" process mpA scattering. Dashed Unes represent 
hard interactions. 
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